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Abstract 

Due to the deflection of light by density fluctuations along the line of sight, weak lensing is 
fNj ■ an unavoidable systematic uncertainty in the use of type la supernovae (SNe la) as cosmological 

lO ' distance indicators. We derive the expected weak lensing signatures of SNe la by convolving the 

m ' 

^O ' intrinsic distribution in SN la peak luminosity with magnification distributions of point sources. 

We analyze current SN la data, and find marginal evidence for weak lensing effects. The statistics 

is poor because of the small number of observed SNe la. Future observational data will allow 

Q-f unambiguous detection of the weak lensing effect of SNe la. The observational signatures of weak 

lensing of SNe la that we have derived provide useful templates with which future data can be 
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PACS numbers: 



I. INTRODUCTION 

The use of type la supernovae (SNe la) as cosmological distance indicators has become 
fundamental in observational cosmology [6, 9, 14, 19-21]. Although SNe la can be calibrated 
to be good standard candles [15, 18], they can be affected by systematic uncertainties. These 
include possible evolution in the intrinsic SN la peak brightness with time [4], weak lensing 
of SNe la [2, 5, 7, 8, 11, 13, 16, 23, 25, 27], and possible extinction by gray dust [1]. 

Weak lensing effect is an unavoidable systematic uncertainty of SNe la as cosmological 
standard candles, simply because there are fluctuations in the matter distribution in our 
universe, and they deflect the light from SNe la (causing either demagnification or magnifi- 
cation). 

In Sec. 2, we derive the expected weak lensing signatures of SNe la by convolving the 
intrinsic distribution in SN la peak luminosity with magnification distributions of point 
sources. In Sec. 3, we use current SN la data to show that weak lensing effect may have 
already begun to set in. Sec. 4 contains a brief summary and discussions. 

II. SIGNATURES OF WEAK LENSING 

The observed flux from a SN la is 

/ = fiL int , (1) 

where L int is the intrinsic brightness of the SN la, and /z is the magnification due to inter- 
vening matter. Note that fi and Li nt are statistically independent. The probability density 
distribution (pdf) of the product of two statistically independent variables can be found 
given the pdf of each variable (for example, see [10]). 
We find that the pdf of the observed flux / is given by 
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P(f)= -r^9(L**)p(T-), (2) 
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where g(L int ) is the pdf of the intrinsic peak brightness of SNe la, p(fi) is the pdf of the 
magnification of SNe la. The upper limit of the integration, L™ x = f /fi m i n , results from 

/* J I ■'-'int — H'min- 

A definitive measurement of g(Li nt ) will require a much greater number of well measured 
SNe la at low z than is available at present. Since g(L int ) is not sufficiently well determined 



at present, we present our results for two different g(Lj ni )'s: Gaussian in flux and Gaussian 
in magnitude. 

The p(/i)'s can be computed numerically using cosmological volume N-body simulations 
(see for example, [2, 16, 24, 25]). We can derive the p(p) for an arbitrary cosmological model 
by using the universal probability distribution function (UPDF) of weak lensing amplification 
[27, 30], with the corrected definition of the minimum convergence [33], 

^n(z) = -- —^- J q dz -^- -^ [X(Z) ~ \{Z )], (3) 

where r(z) is the co moving distance in a smooth universe, 

e{z) = y/n m (i + zf + n x p x (z)/p x (o) + n k (i + z y, 

with px(z) denoting the dark energy density. The affine parameter 

A « - ^ [ ixTwm- (4) 

Note that p m in — 1/(1 — ^mm) 2 - We have used a modified UPDF [34], with the corrected 
minimum convergence and extended to high magnifications. The numerical simulation data 
of p(p) is converted to the modified UPDF, pdf of the reduced convergence 77, 
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where 77 = 1 + (p — l)/\p m in — 1|- The parameters of the UPDF, r] pea k, w, and q are only 
functions of the variance of 77, £ v , which absorbs all the cosmological model dependence. 
The functions r) pea k(C, v ) , w(£, v ), and q(C, v ) are extracted from the numerical simulation data. 
For an arbitrary cosmological model, one can readily compute ^ [30], and then the UPDF 
(and hence p(p)) can be computed. 

Fig.l shows the prediction of the observed flux distributions of SNe la, with magnification 
distribution p(p) given by a ACDM model (f2 m =0.3, ^a = 0.7) at z = 1.4 (top panel) and 
z = 2 (bottom panel) respectively. We have assumed that the intrinsic peak brightness 
distribution, g(L int ), is Gaussian with a rms variance of 0.193 (in units of the mean flux, 
and chosen to be the same as the 0.02 < z < 0.1 subset of the Riess 2004 sample). 

Fig.2 is the same as Fig.l, except here we have assumed that g(Li nt ) is Gaussian in 
magnitude, with a rms variance of 0.213 mag (chosen to be the same as the 0.02 < z < 0.1 
subset of the Riess 2004 sample). 



Clearly, there are two signatures of the weak lensing of SNe la in the observed brightness 
distribution of SNe la. The first signature is the presence of a non-Gaussian tail at the bright 
end, which is due to the high magnification tail of the magnification distribution. The second 
signature is the slight shift of the peak toward the faint end (compared to the pdf of the 
intrinsic SN la peak brightness), which is due to p(fi) peaking at \x < 1 (demagnification) 
because the universe is mostly empty. As the redshift of the observed SNe la increases, the 
non-Gaussian tail at the bright end will grow larger, while the peak will shift further toward 
the faint end (see Figs. 1-2). 

If the distribution of the intrinsic SN la peak brightness is Gaussian in flux, the dominant 
signature of weak lensing is the presence of the high magnification tail in flux. If the 
distribution of the intrinsic SN la peak brightness is Gaussian in magnitude, the dominant 
signature of weak lensing is the shift of the peak of observed magnitude distribution toward 
the faint end due to demagnification. This is as expected, since the magnitude scale stretches 
out the distribution at small flux, and compresses the distribution at large flux. 

III. EVIDENCE OF WEAK LENSING IN CURRENT SUPERNOVA DATA 

We use the Riess sample [20] to explore possible weak lensing in current SN la data, as 
this sample contains the largest number of SNe la at z > 1 that are publicly available. 

Our high z subset consists of 63 SNe la from the Riess sample with with 0.5 < z < 
1.4. Our low z subset consists of 47 SNe la from the Riess sample with 0.02 < z < 0.1. 
Table 1 shows the redshift distribution of the 63 SNe la in the high z subset. To enable 
meaningful comparison of the high z and low z samples, the bestfit cosmological model 
has been subtracted from the brightnesses of all the SNe la in both samples. The bestfit 
cosmological model (found by allowing the dark energy density to be a free function given 
by 4 parameters) has been obtained via flux-averaging of the gold set of 157 SNe la of the 
Riess sample, combined with CMB and galaxy survey data [32], hence it should have very 
weak dependence on the mean brightnesses of the low z and high z SN la samples. 

Table 1 

Redshift distribution of SNe la in the high z subset 



z 


[.5, .7] [.7, .9] [.9, 1.4] 


number 


27 21 15 



Fig.3 shows the flux distributions of 63 SNe la with 0.5 < z < 1.4 (top panel) and 47 SNe 
la with 0.02 < z < 0.1 (bottom panel). The bestfit cosmological model obtained by [32] 
using flux-averaging (which minimizes weak lensing effect, see [29, 31]) has been subtracted. 

Fig.4 is the same as Fig.3, except here we have fitted the distribution of the low z SNe 
la to a Gaussian in magnitude, and have binned the high z SNe la in magnitude as well for 
comparison. 

The top panels of Figs. 3 and 4 show the predicted distributions of SN la peak brightness, 
obtained by convolving the bestfit Gaussian distribution at low z (dotted line) with p(fi) 
for the bestfit cosmological model with z = 1.4 (solid line), and the Poisson noise expected 
from the finite number of SNe la in each bin (dashed lines). 

Clearly, the distribution of the low z SNe la is consistent with Gaussian in both flux and 
magnitude, while the high z SNe la seem to show both signatures of weak lensing: the high 
magnification tail at bright end and the demagnification shift of the peak toward the faint 
end. Within the uncertainties of the Poisson noise, the flux and magnitude distributions of 
the high z SNe la are roughly consistent with the upper bound set by the maximal expected 
amount of weak lensing in the best fit model. 

Both the presence of the bright end tail and the shift of the peak toward the faint end 
are more pronounced than the predictions based on standard lensing magnifications (see 
Figs. 3-4). However, since the number of SNe la at high z is still small, the difference may 
have resulted from statistical fluctuations (Poisson noise). 

Table 2 compares the high z and low z SNe la plotted in Figs. 3 and 4. The mean 
brightness and skewness of the distributions have been calculated for both flux and mag- 
nitude distributions. The rms variance of the skewness S3 for a Gaussian distribution is 
<7<f 3 ~ \/6/N [17], where N is the number of SNe la in the subset. 



Table 2 

Comparison of high z and low z SNe la 



z of SNe 


(IT 


S3 ± erg (flux 


) (mag) b 


S 3 ±ag(mng) 


.5 <z < 1.4 


1.115 


2.12 ±0.31 


-.0763 


-1.06 ±0.31 


.5 < z < 1.4 C 


1.055 


0.66 ±0.32 


-.0352 


-0.25 ±0.32 


.02 < z < .1 


1.096 


0.60 ±0.36 


-.0805 


-0.19 ±0.36 



a in units of the flux in the bestfit cosmological model obtained via flux-averaging in [32] . 



6 defined to be — 2.5 \og(flux), where flux is the dimensionless flux defined above. 
c excluding the three brightest SNe listed in Table 3. 

The high z and low z SNe la differ little in mean brightness; this is as expected since the 
bestfit cosmological model obtained using the gold set of 157 SNe la of the Riess sample 
(together with CMB and galaxy clustering data) has been subtracted. 1 In the absence 
of weak lensing, one would expect that the low z and high z SN la samples have similar 
distributions, hence similar values of the skewness S3. If the intrinsic distribution of SN la 
brightnesses were Gaussian in flux (magnitude), one would expect S3 = ± <7<j for the data 
in flux (magnitude). Table 2 shows that: (1) The high z SN la sample has a significantly 
larger skewness S3 (for both flux and magnitude distributions) than the low z sample; (2) 
The large skewness of the high z sample is primarily due to the three brightest SNe la. Both 
these points are consistent with the signatures of weak lensing discussed in Sec. 2. 

Table 3 lists the three brightest SNe la in the bright end tails of Figs. 3 and 4. All three 
SNe la are in the "gold" sample of [20] . The last column in Table 3 lists the possible range of 
magnification ji for each SN la, [(/ — df)/(f), (f + df)/{f)], with df/f = sinh (<7 Mo In 10/2.5), 
and (/) given by the low z subset. Note that the possible magnification of SN1998I has a 
very large uncertainty, because its distance modulus has a very large uncertainty in the Riess 
sample: \x§ = 42.91 ± 0.81, which correspond to a flux uncertainty of about 80%. 

Table 3 
Three brightest SNe la 



SN 


z no 


<Vo 


possible ji 


SN1997as 0.508 41.64 0.35 


[1.42, 


2.78] 


SN2000eg 


0.540 41.96 0.41 


[1.10, 


2.50] 


SN1998I 


0.886 42.91 


0.81 


[0.44, 


4.40] 



Note that the bestfit cosmological model [32] was obtained using the gold set of 157 SNe 
la of the Riess sample [20], flux- averaged and combined with CMB and galaxy clustering 
data. The mean brightnesses of the low z and high z samples considered in this paper were 
not used in deriving the bestfit cosmological model (which is given by 6 parameters). The 
fact that the low z and high z samples have about the same mean brightnesses shows the 



1 Flux-averaging has been used in obtaining the bestfit model in order to minimize the bias due to weak 
lensing [32]. 
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validity of the bestfit cosmological model, since weak lensing does not change the mean 
brightness of teh high z sample due to flux conservation. 

The second row in Table 2 shows that if we exclude the three brightest SNe listed in 
Table 3, the skewness of the high z sample drops to about the same as that of the low z 
sample. However, the mean brightness of the high z sample drops below that of the low z 
sample, such that the high z sample is about 4% fainter in flux than the low z sample. This 
suggests that the three brightest SNe in the high z sample are probably not outlyers, but 
may belong to the high magnification tail of p(/j,). 2 

Finally, we do a Kolmogorov-Smirnov test to assess whether the low z and high z SNe la 
could be from the same brightness distribution (i.e., no weak lensing of the high z sample). 
Table 4 shows the Kolmogorov-Smirnov test of the low z versus high z SN la sample, 
and the high z SN la sample compared to a Gaussian in magnitude (with a mag = 0.213) 
convolved with p(fi) at z = 1.4 (as plotted in Fig. 4a). We have chosen the later distribution 
for comparison with the high z sample, because the low z sample appears slightly more 
Gaussian in magnitude (with smaller S3, see Table 2). 

Table 4 
Kolmogorov-Smirnov test 



subtracted cosmological model 


bestfit model 


V^ ^m _ 


= 0.27, n A = 


= 0.73) 




D prob. 




D 


prob. 


low z and high z sample 


0.153 0.522 




0.184 


0.291 


high z sample and p(mag) a 


0.096 0.585 




0.104 


0.482 



a a Gaussian in magnitude (with a ma g = 0.213) convolved with p(/i) at z — 1.4, see Eq.(2). 

The Kolmogorov-Smirnov test gives D, the maximum value of the absolute difference 
between two cumulative distribution functions, and prob, the probability that D > observed. 
Small values of prob show that the two distributions are significantly different. Table 4 shows 
that current data do not yield results of high statistical significance, however, the high z SN 



2 Note that weak lensing leads to a continuous probability distribution function (pdf) in magnification with 
a tail at high magnifications which are often associated with strong lensing. However, strong lensing 
usually refers to lensing that results in multiple images of a point source [22] . The high magnification tail 
of a weak lensing pdf corresponds to the bending of light of a point source due to galaxies that results 
in a single magnified image [26], hence it is technically still weak lensing magnification, and not strong 
lensing. 



la sample is more consistent with a lensed distribution than with the low z SN la sample. 
Changing the subtracted cosmological model from the bestfit model to a popular model with 
Q m = 0.27, Q\ = 0.73 does not change our results qualitatively. 

IV. SUMMARY AND DISCUSSION 

We have derived the expected weak lensing signatures of SNe la in the distribution of 
observed SN la peak brightness, the presence of a high magnification tail at the bright end 
of the distribution, and the demagnification shift of the peak of the distribution toward the 
faint end (see Figs. 1-2). 

Our method is complementary to that of [12, 35]; they use the correlation between fore- 
ground galaxies and supernova brightnesses to detect weak lensing, while our method only 
uses the statistics of supernova brightnesses and does not depend on the observation of 
foreground galaxies. 

We have compared 63 high z SNe la (0.5 < z < 1.4) with 47 low z SNe la (0.02 < z < 0.1) 
from the Riess sample [20]. We find that the observed flux and magnitude distributions of 
the high z sample are roughly consisitent with the maximal expected amount of weak lensing 
magnification (see Figs. 3-4), within the Poisson noise due to the small number of SNe la in 
each bin. 

We have identified the three brightest SNe la in the high z subset (0.5 < z < 1.4) of 
the Riess sample, and estimated a possible range of magnification for each SN la (Table 3). 
Observational follow-up of the regions near these SNe la may show whether these SNe la 
have indeed been magnified. Note that selection effects should not be important here, since 
the three brightest SNe are at intermediate redshifts (where the fainter SNe are not close to 
the detection limits of the surveys). 

Our results are consistent with those of [3, 35]. [35] found that brighter SNe la are 
preferentially found behind regions which are overdense in foreground galaxies, as expected 
in weak lensing. [3] found tentative evidence for a deviation from the reciprocity relation 
between the angular diameter distance and the luminosity distance, which could be due to 
the brightening of SNe la due to lensing. 

Although the observed flux/magnitude distribution of the high z SN la sample deviates 
from a Gaussian in ways that are qualitatively consisitent with the weak lensing effects (al- 



8 



lowing for Poisson noise), it is possible that these deviations could be due to other unidenti- 
fied systematic effects. However, it is important to note that if we remove the three brightest 
SNe la in the high z sample, the distributions become more Gaussian, but the mean becomes 
biased (see Table 2). Therefore, these three SNe la might not be outliers, and should be 
included in the data analysis. 

Fortunately, the non-Gaussianity of the high z SN la flux/magnitude distribution (re- 
gardless of its origin - weak lensing or some other systematic effect) does not seem to alter 
the mean of the distribution (compared to the low z sample). This means that we can, and 
should, use flux averaging in order to obtain unbiased estimates of cosmological parameters 
[29, 31]. 

Our results explain the difference between the cosmological constraints found by [20] and 
[32] for the same model assumptions (see Fig. 10 of [20] and Fig. 2(a) of [32]). [32] used 
flux- averaging in their likelihood analysis; [20] did not. 

As more SNe la are discovered at high z, it becomes increasingly important to minimize 
the effect of weak lensing (or other non-Gaussian systematic effects that conserve flux) by 
flux-averaging [29, 31] 3 in using SNe la to probe cosmology. 

It seems that weak lensing effects, or some other systematic effect that mimics weak 
lensing qualitatively, may have begun to set in (see Figs. 3-4 and Tables 2 & 4). However, 
the statistics is poor because of the small number of observed SNe la. Future observational 
data from current and planned SN la surveys will allow unambiguous detection of the weak 
lensing effect of SNe la. The observational signatures of weak lensing of SNe la that we have 
derived (see Figs. 1-2) provide useful templates with which future data can be compared. 

Public software: A Fortran code that uses flux-averaging statistics to compute the 
likelihood of an arbitrary dark energy model (given the SN la data from [20]) can be found 
at http : / /www. nhn.ou.edu/ ~ wang/SNcode/. 



If the distribution of the intrinsic peak brightness of SNe la is Gaussian in magnitude, then flux-averaging 
would introduce a small bias of — o^ In 10/5[28], which needs to taken into account in the data analysis. 
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FIG. 1: Prediction of the observed flux distributions of SNe la for magnification distribution p(fi) given by 
a ACDM model (J7 m =0.3, £l\ = 0.7) at z = 1.4 (top panel) and z — 2 (bottom panel) respectively. We have 
assumed that the intrinsic peak brightness distribution, g(Li nt ), is Gaussian with a rms variance of 0.193 
(in units of the mean flux) . 



12 



1 - 
p(mag) ; 

0.5 - 



p(mag) ; 

0.5 - 







"i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — r 
Gaussian with a m . =.213 

mag 



p(/x) [ACDM, n m = .3, z=1.4] 

Gaussian convolved with p(/x) 



o hi--!- i--{-r-:-i-i- 
l 




Gaussian with a m „ =.213 

mag 



P(aO [ACDM, n m = .3, z=2] 

Gaussian convolved with p(/x) 



J™L L__U_I I l_ 




-1.5 -1 



■0.5 

Amag 



0.5 



J L 



1 



FIG. 2: Same as Fig.l, except here we have assumed that the intrinsic peak brightness distribution, g(Li n t), 
is Gaussian in magnitude, with a rms variance of 0.213 mag. 
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FIG. 3: The flux distributions of 63 SNe la with 0.5 < z < 1.4 (top panel) and 47 SNe la with 0.02 < z < 0.1 
(bottom panel). The bestfit model obtained by [32] using flux- averaging (which minimizes weak lensing 
effect) has been subtracted. Clearly, the distribution of the low z SNe la is consistent with Gaussian, while 
the high z SNe fa seem to show both signatures of weak lensing (high magnification tail and demagnification 
shift of the peak to smaller flux). The error bars indicate the Poisson noise of the weak lensing prediction 
(solid curve). 
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FIG. 4: Same as Fig. 3, except here we have fitted the distribution of the low z SNe la to a Gaussian in 
magnitude, and have binned the high z SNe la in magnitude as well for comparison. 
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